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We have studied the neutral Higgs boson (NHB) contributions to the pure penguin process B — > 
4>Ks in MSSM with middle and large tan/3 (say, > 8). We show that the ct B order hadronic 
matrix elements of NHB induced operators can make sizable effects on both the branch ratio and 
time dependent CP asymmetry S^k- Under the all relevant experimental constraints, the Higgs 
mediated contributions to S^k alone can provide a significant deviation from SM and, in particular, 
lead to a negative S^k which is reported by BaBar and Belle. 



Among charmless hadronic B decays, the rare decay B —> <&Ks is one of channels which provide a powerful testing 



o 
o 



ground for new physics. The reason is simply because it has no tree level contributions in the standard model (SM). 
The decay has been studied in SM and the predicted branching ratio (Br) in the perturbative QCD framework is 
(3.5-4.3) x 10" 6 Q by using the BBNS approach @,||] or about 10 x 1CT 6 by using the H.-n. Li et al's approach 0, 
which is somewhat smaller than or consistent with the measured value (8.4^ i) x 10 -6 , the current world average 
from BaBar 0] and Belle 0. 

The recently reported measurements of time dependent CP asymmetries in B — > $>Ks decays by BaBar 6] and 
Belle lead to the error weighted average 
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with errors added in quadrature*. 

Here the time-dependent CP-asymmetry S^k is defined by 



<N ; S^ K = sin(2/3($Ks)W = -0.39 ± 0.41 (1) 
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q ; a^xit) = -C^k cos(AM B oi) + S^k sin(AM B oi), (2) 

C^) ' where 

CD , with A^/s being 



,q\ A(B° -> <j>K s ) 
4>k = - 7750 . TtFT\ ■ ( 4 ) 
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In the SM the above asymmetry is related to that in B — > J/^S/Ks by 

sin(2/3($Jf)) = sin(2f3(J/^K)) + 0(X 2 ) (5) 

where A ~ 0.2 appears in Wolfenstein's parameterization of the CKM matrix and 

sm{2(3(J/^K s ,L))worid-ave = 0.734 ± 0.054. (6) 

Therefore, JIJ violates the SM at the 2.7 a deviation. 

Obviously, the impact of these experimental results on the validity of CKM and SM is currently statistics limited. 

However, they have attracted much interest in searching for new physics mm m. 

In the letter we inquire the possibility to explain the 2.7 a deviation in R-parity conservative supersymmetric models, 
e.g., the minimal supersymmetric standard model (MSSM). There are mainly two new contributions arising from the 
QCD and chromomagnetic penguins and neutral Higgs boson (NHB) penguins with the gluino and squar k propa gated 
in the loop. The QCD and chromomagnetic penguin contributions to b — > sss have been analyzed in Refs.|l2Ul3|. The 
NHB penguin contributions have been discussed in Ref.[13j. However, the conclusion on NHB contributions, S$k can 



* 2003 data is -0.96 ± 0.50+°;^ by Belle || and it is +0.45 ± 0.43 ± 0.07 by BaBar 0. 
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not smaller than 0.71, in Ref.[l3j is valid only in some special cases, e.g., the extended minimal flavor violation (MFV) 
scenarios 13J with the naive factorization of hadronic matrix elements of operators and the decoupling limit in MSSM, 
and it is not valid in (general) MSSM, as we shall show. We concentrate on the Higgs penguin contributions in MSSM 
in the letter. An interesting fact is that the Higgs penguins contribute to B — > (j>Ks but not B — > J/tpKs at the a s 
order of hadronic matrix element calculations due to the non-match of color and flavor. At the tree level of hadronic 
matrix element calculations (i.e., the naive factorization) their contributions to the branch ratio of B — > J/ipKs (the 
time dependent CP asymmetry Sj/^k) are very small (negligible). Therefore, we do not need to worry about the 
effects of the Higgs penguins on Sj/^r. It is shown that the Higgs mediated contributions in the case of middle and 
large tan/3 (say, > 8) are important under the constraints from the experimental bounds of B s — > /i + [J,~ and AM S 
as well as b — > s 7, s g, B — ► X s u, + n~ and consequently can have significant effects on B — > §>Ks due to the tan 2 (3 
enhancement of Wilson coefficients of the NHB induced operators in MSSM [Til ITHj l . Our results show that the Higgs 
mediated contributions to Sqk alone can provide a significant deviation from the SM and a possible explanation of 
the 2.7 a deviation, satisfying all the relevant experimental constraints. 

The effective Hamiltonian induced by neutral Higgs bosons can be written as 

H n e t = %{-h) £ (C l Q l + C' l Q' 1 ) + h.c, (7) 

v Z z=ll,. ...16 



where Qn to Qm, the neutral Higgs penguins operators, are given by 

Qll(13) 



{sb) S+P J2 -%2)s-(+)P: 
q mb 



?12(14) = (sibjjs+p^ — (qj qj)s-(+)P , 

TTlfj 



Q 15 = sa^(l + 75 )bJ2^q a^(l+ l5 )q, 
q 

Qi6 = s t a" v (l + 75) bj V ^-qj a„ u (l + 75) $ , (8) 

q 

where {qiq2)s±p — 5i(l i 7s)?2 ■ The operators Q^s are obtained from the unprimed operators QiS by exchanging 
/. • • //. 

There are new sources of flavor violation in MSSM. Besides the CKM matrix, the 6x6 squark mass matrices are 
generally not diagonal in flavor (generation) indices in the super-CKM basis in which superfields are rotated in such a 
way that the mass matrices of the quark field components of the superfields are diagonal. This rotation non-alignment 
in the quark and squark sectors can induce large flavor off-diagonal couplings such as the coupling of gluino to the 
quark and squark which belong to different generations. These couplings can be complex and consequently can induce 
CP violation in flavor changing neutral currents (FCNC). It is well-known that the effects of the primed counterparts of 
usual operators are suppressed by m s /mb and consequently negligible in SM because they have the opposite chiralities. 
However, in MSSM their effects can be significant, since the flavor non-diagonal squark mass matrix elements are free 
parameters. We have calculated Wilson coefficients of Qi and Q\ coming from the gluino-squark loop contributions 
at the m w scale and their explicit expressions are given in R ef . II 51 . We evolve the Wilson coefficients to the low scale 
(i ~ nib using the one loop QCD anomalous dimensions 0, Il7lll8j . 

The effective Hamiltonian results in the following matrix element for the decay B — > <f>K s 

(K s cj>\HT\B) = % E A f { K ^W eU \ B ) > ( 9 ) 

* p—u,c 

where T p neu is given by 

j-neu =a n-(s6) v _ A ®(s S ) v _ A 



t Strictly speaking, the sum over q in expressions of Qi should be separated into two parts: one is for q=u, c, i.e., upper type quarks, the 
other for q=d, s, b, i.e., down type quarks, because the couplings of upper type quarks to NHBs are different from those of down type 
quarks. In the case of large tan j3 the former is suppressed by tan -1 j3 with respect to the latter and consequently can be neglected. 
Hereafter we use, e.g., Cjj to denote the Wilson coefficient of the operator Qn = (sb)g+p ^-(cc)g_p. 
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In Eq.dJ 
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where 



r 1 /7 T r 1 
G% = I — $ (x), GF{a,x)= / dtln[a-a:«] 
Jo x Jo 



GFJs) = [ dx^^GF( 
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(13) 



with x = 1 — x and = 6xx in the asymptotic limit of the leading-twist distribution amplitude. In calculations 

we have set m U: d = and neglected the terms which are proportional to m^/ffij in Ea. (|12fl . We have included only 
the leading twist contributions in Eq. l|ll|) . In Eci. H10[) a' is obtained from a, by substituting the Wilson coefficients 
CjS for CjS. Here we have used the BBNS approach Q to calculate the hadronic matrix elements of operators. 
Notations not explicitly defined are the same as those in Ref.Q. 

Before numerical calculations a remark is in place. From Eq. 112fl we can see that the large contributions to 

matrix elements of the operators Q\'\ i=ll,...,16 arise from penguin contractions with b quark in the loop. It is the 
contributions which make the effects of Higgs mediated mechanism sizable at the a s order. At the fj, = m w scale only 
non zero Wilson coefficients are G^ 13 , and C[§ 16 (/Lt) are obtained from Cj2(rn w ) due to the operator mixing under 
renormalization. Moreover, because of the mixing of Qi, i=13,...,16, onto Q77.sc/, Csg(^) (m ~ m b) can be significantly 
enhanced [Tsj] . which can have large effects on both Br and S^k, as can be seen from the SM decay amplitude of 
B -> <f>K s . 

In numerical calculations we use the SM decay amplitude of B — > 4>K S given in Ref.Q. We calculate the evaluation 
of Wilson coefficients to the next-to-leading for the SM operators and to the leading order for <5n,...,i6 and their 
mixing with the SM operators. In order to see how large the effects of NHBs can be we switch off all other SUSY 
contributions except for those coming from NHB mediated penguin operators. We show model-independently in FIG. 
1 Sqk versus the phase of C^(rn w ) with fixed values of \C\^(rn w )\, setting Cn(m w ) and C' lx 13 (m w )=0 in order to see 
the essential feature of the a s order corrections due to the Higgs mediated mechanism. We find that S^k can reach 



FIG. 1: S^k versus the phase of C\3(m w ) for fixed |Ci3(m„)| for (a) fi = mi, (b) /i = rrib/2. The solid, dotted and dashed 
curves correspond to |Ci3(m m )|= 0.05, 0.15 and 0.25 respectively. 
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FIG. 2: The correlation between S^k and Br(_B — » </>A" s ) for including only SM and NHB contributions, (a) is the tree level 
result for /i = mi,, (b) is including the a s corrections. In Fig. 2b green and red dots denote the results for /i = mb and 
fi = mt/2 respectively. Current la bounds are shown by the dashed lines. 

a minus value if |Ci3(m u) )| is equal to or larger than 0.10 for /i = m& and 0.11 for [i — mi,/2. S^k can be negative 
in quite a large range of the phase of C13 for a larger value of \Ci^{m w )\. S^k is dependent of but not sensitive to 
the characteristic scale \i of the process. For a specific model, e.g., MSSM which in the letter we are interested in, we 
calculate all relevant Wilson coefficients c['i 13 {m w ) in reasonable regions of parameter space with middle and large 
tan/5 and present numerical results in FIG. 2 below. 

We impose two important constraints from B — > X s "f and B s — > Considering the theoretical uncertainties, 

we take 2.0 x 10~ 4 < Bi(B — > X s -f) < 4.5 x 10~j, as generally analyzed in literatures. The current experimental 
upper bound of Br(i3 s — * is 2.6 x 10 -6 [l!j. Because the bound constrains \Cq 1 — Cq.| (i=l, 2)* we can 

have values of |CqJ and \C'q. \ larger than those in constrained MSSM (CMSSM) with universal boundary conditions 
at the high scale and scenarios of the extended minimal flavor violation in MSSM in which \C'q.\ is much smaller 
than |CqJ. We require that predicted Br of B — > X s ^ + [i~ falls within 1 a experimental bounds. We also impose 
the current experimental lower bound AM S > 14.4p.s~ 1 and experimental upper bound Br(_B — > X s g) < 9%. In 
numerical analysis we fix rrig — rriq — 400GeV and tan/3 = 30. We vary the NHB masses in the ranges of 91GeV < 
m h < 135GeV,91GeV < m H < 200GeV with m h < m H and 200GeV < m A < 250GeV for the fixed mixing angle 
a = 0.6, 7r/2 of the CP even NHBs and scan 8^ A in the range \o2% A \ < 0.1 which arises from the constraints of b — > sj 
and AM S (A = L, R). Here is the parameter in the usual mass insertion approximation [21.] and its definition 

can be found in Refs. [15L l2l| . We perform calculations of hadronic matrix elements of the Higgs penguin operators 

at the a® order (tree level), i.e, in the naive factorization formalism, and to the a s order. The results for tree 
level (i.e., the naive factorization) and to the a s order are shown in FIGs. 2a and 2b, respectively. Fig. 2a is given 
for (i — TO;,. In Fig. 2b green and red dots denote the results for jj, = mj, and /j, = m&/2 respectively. We find that 
although S$k can have a sizable deviation from SM but it can not be smaller than 0.6 if one use the naive factorization 
to calculate hadronic matrix elements. However, S$k for both [i — m\, and fi — mf,/2 can have a significant deviation 
from SM when we include the a s corrections. We find that there are regions of parameters where S^k fells in la 
experimental bounds and Br is smaller than 1.6 x 10~ 5 . The FIGs. 2a and 2b are plotted for a = ir/2. The similar 
figures are obtained for a = 0.6. We stress that at tree level S^k can have a sizable deviation from SM because we 
can have values of \Cq < \ and \Cq. | larger than those in CMSSM with universal boundary conditions at the high scale 
and the extended MFV scenarios in MSSM by a factor of 3 or so, as pointed out above. It is instructive to note 
that because the LR and RL mass insertions do not contribute to the Higgs penguin operators, the Higgs mediated 



i C*q 2 are the Wilson coefficients of the operators Qi \ which are Higgs penguin induced in leptonic and semileptonic B decays and their 

definition can be found in Ref. L'll . By substituting the quark-Higgs vertex for the lepton-Higgs vertex it is straightforward to obtain 
Wilson coefficients relevant to hadronic B decays. 
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FIG. 3: The correlation between S$k an d Br(B — > </>iG) for (a) including only SM and NHB contributions and only a LL 
insertion, (b) all contributions in MSSM. Green and red dots denote the results for /j, = nib and fi = mt/2 respectively. Current 
la bounds are shown by the dashed lines. 



mechanism probes the LL and RR insertions, in contrast with the QCD and chromomagnetic penguins whose effects 
are too small to alter S^k significantly in the case of only LL and RR insertions |l3| ^ • We have also carried through 
the calculations in the middle tan/3 case, tan/3 = 8, within the range l^^ 1 ! — 1 an d the results are similar to those 
in the large tan/3 case. 

Next we consider the case of large m,A- In the limit M\ M§ (for Ma > 300GeV), the charged, heavy CP-even 
and CP-odd neutral Higgs bosons are nearly mass degenerate, M H ± ~ Mh — Ma, sin(/3 — a) approaches 1, and the 
properties of the lightest CP-even neutral Higgs boson h are almost identical to those of the SM Higgs boson (so called 
the decoupling limit). In this case, C'['^ approaches to zero if rribAb <C m~. So S$k for /i ~ nib can not be smaller than 
0.7 under the constraint from B s — > ji + [i~ . However, as shown in Ref. |22j |. the decoupling limit can be relaxed in 
some regions with M\ ^S> M§ of the parameter space which are allowed by experiments of (g — 2)^, b — + sj and lower 
bounds of Higgs and sparticle masses, due to the large off-diagonal scalar top and scalar bottom mass matrix elements 
contributing to the Higgs sector by radiative corrections. In the regions, the charged, heavy CP-even and CP-odd 
Higgs bosons are not mass degenerate, and sin 2 (/3 — a) can be damped from 1. We take two typical set of parameters 
in the regions given in the TABLE 1 Case B of Ref.[I2 to calculate Wilson coefficients c\'\m w ), i = 11,13,77. 
Corresponding numerical results are that S^k = —0.53 and —0.19, Br = [0.77 and 0.98] x 10~ 5 , respectively. 

We now consider the case of only a LL insertion. Then C'n 16 = and consequently the experimental bound of 
B s — > /i + fj,~ constrains |CqJ (i=l, 2). We scan (Jfij in the range | ^23^ \ < 0.1 with other parameters same as above 
under all relevant constraints. The result is shown in Fig. 3a. We find that S^k can significantly deviate the SM 
value but it can not be smaller than 0.4. As pointed out above, if using the naive factorization of hadronic matrix 
elements, the minimal value of S$k is 0.7, in contrast with that of including the a s corrections of hadronic matrix 
elements. The similar result is obtained for the case of only a RR insertion. 

Finally we switch on other SUSY contributions of which the contributions of the chromomagnetic penguin with 
gluino-down type squark in the loop are dominant. We scan \S2^ L ' RR \ < 0.1, with values of other parameters the same 
as those above under the constraints from b — > sj and AM S as well as other relevant experiments mentioned above. 
The result including both the chromomagnetic and QCD penguin contributions and NHB penguin contributions is 
shown in FIG. 3b. One can see from the figure that there are regions of parameters where S^k falls in la experimental 
bounds and Br is smaller than 1.6 x 10 -5 . Comparing FIG. 3b with FIG. 2b, one can see that the regions of parameter 
space in the case including both contributions are larger than those when including only the NHB contributions. 

Our results show that both the Br and S$k of B — > <f)Ks are dependent of the scale fi , which can be understood 
as follows. Let us look at, e.g., the contributions coming from <3i3,...,i6 in a™ e " (see, Eq. I|12|l '). we have: 

dal eu 



dln/i 



% & (C 13 2 In » - 4C 15 2 In M - 8C 16 2 In M )G° 



dln/i 

4tt AN r 



[(C13 2 - 4C 15 2 - 8C 16 2))]G° + a* hif! terms + (14) 



' It is possible that effects of the QCD and chromomagnetic penguins are significant to alert S^k greatly for only a LL or RR insertion 
because there can be an induced LR or RL insertion if /i tan j3 is large enough ll 'I . 
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where denotes the terms coming from Qn and Q' i} i=ll,13,...,16. The first term in Eq. I|14fl is cancelled by 

rf(-2C 8g G°) 
din ji 

= £4^[- 2(Ci3 - 4C,i5 - 8Ci6) ]^ + -' 

where — 2Cg ff G° is the contribution coming from the chromomagnetic dipole operator and denotes the terms 
independent of Cia,...,^. However, the a^ln/i terms in Eq. (|T4"|l is left. Although their appearance is of a 2 s order, 
their magnitude indeed are of a s order. Similar analysis can be down for the other terms in a\ eu . 

It is necessary to make a theoretical prediction in SM as precision as we can in order to give a firm ground for 
finding new physics. For the purpose, we calculate the twist-3 and weak annihilation contributions in SM using the 
method in Ref. |23| by which there is no any phenomenological parameter introduced^. The numerical results show 
that the annihilation contributions to Br(B — > 4>K) are of order 10~ 8 which are negligible, the twist-3 contributions 
to Br are also very small, smaller than one percent, and both the annihilation and twist-3 contributions to S^k are 
negligible. 

In summary, we have shown that the NHB contributions in MSSM can have a significant effect on both the Br and 
Sa,k of B — > 4>Kg. We can have large Wilson coefficients (comparable with those in SM) of the NHB penguin induced 

operators Ql'\ in some regions of the parameter space in MSSM under all relevant experimental constraints, which 
leads to that even at the tree level of hadronic matrix element calculations the time dependent CP asymmetry S^k 
can deviate from SM sizably. Due to the large contributions to the hadronic elements of the operators at the a s order 
arising from penguin contractions with b quark in the loop, both the Br and S^k are sizably different from those in 
SM. The mixing Qi3,...,i6 with Qg g makes the Wilson coefficient Cs g (/i) (/i ~ to&) enhanced a lot, which has large 
effects on both S^k and the Br. In some regions of the parameter space in MSSM, the NHB contributions alone can 
lead to negative S^k which is reported by BaBar and Belle in 2002 and even make S^k agreed with 2003 Belle data 
in a la experimental bounds. Our results show that both the Br and S^k of B — * 4>K$ are dependent of the scale fi. 
The sizable scale dependence, which implies hadronic uncertainties up to the a s order, comes mainly from the 0(a s ) 
corrections of hadronic matrix elements and also from leading order Wilson coefficients C- , i= 11,... ,16. However, 
despite there are hadronic uncertainties, the conclusion that the NHB contributions in MSSM can have a significant 
effect on both the Br and S$k of B — ► <f>Ks can still be drawn definitely. 
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